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Abstract 

We show that the open string worldsheet description of brane decay (dis- 
cussing a specific example of a rolling tachyon background) can be related to a 
sequence of points of thermodynamic equilibrium of a grand canonical ensemble 
of point charges on a circle, the Dyson gas. Subsequent instants of time are 
related to neighboring values of the chemical potential or the average particle 
number N. The free energy of the system decreases in the direction of larger N 
or later times, thus defining a thermodynamic arrow of time. Time evolution 
equations are mapped to differential equations relating thermal expectation 
values of certain observables at different points of thermal equilibrium. This 
suggests some lessons concerning emergence of time from an underlying micro- 
scopic structure in which the concept of time is absent. 
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1 Introduction 

There are a number of examples in string theory in which several dimensions of space 
emerge as an effective description of a system that does not contain them [1-3]. All 
these examples arise at their root from the duality between open and closed strings, 
by equivalently describing the dynamics of a system of D-branes in terms of the open 
string degrees of freedom that quantize them, or in terms of the closed strings that 
quantize the spacetime that the D-branes create. There are no similar examples of 
the emergence of time^- A candidate setting for investigating such an emergence is 
provided by the unstable branes of bosonic and superstring theory [7]. In this case, 
open strings are localized in time [8], and one might ask how the open strings contrive 
to describe closed strings and spacetime at times far from the region of time where 
they are confined. 

Previous work showed that the worldsheet correlation functions of strings in such 
universes can be computed at tree level from an ensemble of U(N) matrices [9-14]. 
As the brane decays to a gas of closed strings, the rank of the matrices contributing to 
the ensemble increases to infinity. In examples of the emergence of space [1-3], closed 

1 For recent discussions of emergence of spacetime and string theory, see [4,5]. For a recent 
discussion from a condensed matter physics perspective, see [6]. 
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string backgrounds emerge from the dynamics of U(N) matrices in the large N limit. 
This suggests that for decaying branes the geometry at late times is "emerging" from 
the dynamics of the large rank matrices appearing in the open string description. 

Here we take a related perspective, by showing that all worldsheet correlation 
functions in decaying brane universes can also be computed at tree level in string 
theory from a simple statistical mechanical model - test charges interacting with a 
Dyson gas, consisting of point charges confined to a circle. Physical time maps into 
(an analytic continuation of) the fugacity in the gas. Thus, by a Legendre transform, 
time becomes related to the mean number of particles (N) in the statistical system. 
(In turn, this is related to the rank of the matrices appearing in the matrix model 
description of the system [9, 12].) The free energy of the Dyson gas decreases with 
N, giving rise to a thermodynamic arrow for time. 

In classical physics, we are used to specifying the state of a system by specifying 
its initial conditions and then allowing the equation of motion to evolve the configu- 
ration through time. Since the Dyson gas arises here from a rewriting of worldsheet 
perturbation theory, the spacetime equations of motion can be derived by requiring 
conformal invariance (perhaps implying criticality) of the statistical system. The 
time evolution equations of spacetime fields then translate into equations relating the 
values of certain statistical observables at different points of thermodynamic equilib- 
rium. 

As a rewriting of worldsheet perturbation theory in a different, statistical language, 
our analysis does not describe a holographic duality in the sense of the AdS/CFT 
correspondence. The latter is a relationship between the spacetime theory of open 
strings (i.e. a quantum field theory) and a spacetime theory of closed strings. In 
principle, we are simply describing a dictionary between two equivalent worldsheet 
formulations. Nevertheless, the appearance of a thermodynamic arrow of time from 
our statistical formulation of the decaying brane is suggestive of some kind of emergent 
phenomenon. The discussion in Sec. 4, attempts to draw lessons about what it could 
mean for time to be "emergent" . 

2 The half S-brane as a Dyson Gas 

2.1 The half S-brane as an ensemble of matrix models 

Bosonic string theory has a spectrum of unstable D-branes that decay by tachyon 
condensation. Sen [7] showed that their homogeneous decay can be described in 
terms of the open string theory on a D-brane by the exactly marginal boundary 
deformation 
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The parameter A has no meaning since it can be absorbed by redefining the origin 
of the target space time coordinate X°. The resulting background is called the half 
S-brane since it represents a spacelike D-brane at t — > — oo that decays into closed 
strings as time passes. 

Worldsheet correlation functions in the background of the decaying brane take the 
form 

Ai= f DX° DX 1 --- DX 25 e~ s fj V a (z a , z a ) (2) 

0=1 

where the action S includes the boundary deformation §T§ and the V a are vertex oper- 
ators constructed from {X°, X} and written as functions of holomorphic coordinates 
on the worldsheet. String scattering amplitudes are obtained by integrating these 
correlation functions over all vertex operator positions (after fixing the location of 
one operator by using the conformal symmetry of the worldsheet). 

Taking V a to be tachyon vertex operators and leaving out the trivial spatial part 
of the computation, the X° dependent piece of the correlator is 

A t ~ f ^e- 5 ^^^"'. (3) 

J o=l 

At non-zero momentum it is possible to choose a gauge in which general on-shell 
closed string vertex operators with finite energy have the form 

V = e ik ° x °V sp , (4) 

where V sp is constructed entirely out of the 25 spatial fields [15]. In this gauge, and 
given that the boundary perturbation only depends on X°, the non-trivial part of 
finite momentum string scattering amplitudes can be written as ([3]). 

The worldsheet correlation function ([3]) can be evaluated on the disk by isolating 
the zero mode x° from the fluctuations as X° = x°+X'°, and expanding the boundary 
perturbation e~ 6Sbdry in a power series. The result is 

J N=0 ' J i=l ^ o=l 

(5) 

Now using the standard Neumann correlator of X'° (see, e.g., Sec. 2.2 and 3.1 of [12]) 
the disk expectation value in (jSj) can be evaluated as 



A l = J dx°e ix °^ a Mx°) (6) 



a<b ab 



where the initial factors involving z a arise from contractions between the vertex oper- 
ators, the integral over x° is the Fourier transformation to total energy, and F isolates 
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all the interesting pieces in the correlator arising from the interaction between the 
vertex operator and the boundary perturbation (fTT) . Ai(x°) is the amplitude written 
as a function of target space time. We find 



Af=0 ' ' i=l i<j io 



't/(iV) 

The second equality arises by recognizing the measure for integration over U(N). 

In this way, worldsheet correlation functions for a decaying brane are expressed in 
terms of an ensemble of U(N) matrices with varying rank N. As a simple example, 
the disk partition function is 

Zo P en = (e- SS ^) = £ } / dU (9) 

(10) 



1 + 2n\e x ° ' 

At early times x° — > — 00 the first few terms in the series ([9]) dominate and hence 
only low rank U(N) matrices contribute. However, as x° increases, and brane decay 
progresses further, matrices of larger rank become progressively more important. The 
series converges when |27rAe x °| < 1. However, the summed expression (jTOl) can be 
analytically continued to future infinity x° — >• 00. The convergence radius x° ~ 
— ln(27rA) might be interpreted as the moment after which the tachyon has completely 
condensed, so that the open string worldsheet is no longer really meaningful. After 
this point, a closed string description is better. Indeed, the function 

is the coefficient of the lowest, closed string tachyon component of the boundary state 
of the decaying brane [7, 9jE 

As another example of an interesting target space observable, consider the energy- 
momentum tensor as a function of target space time x° [9, 18]. It can be decomposed 

as 

T MN (x°) = K[rj MN B(x°) + A MN (x )} , (12) 
where B(x°) = f(x°) is the open string disk partition function, whereas 

A MN (x°) = 2(: dX M (0)8X N (0) :h 5 _ 6mne (13) 



2 For additional discussion and early references, see [11,16,17]. 
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is the disk one-point function of the closed string vertex operator 



V MN {z,z; k) =: OX 



{z)dX N (z)e 



ik-X(z.z) . 



(14) 



at zero momentum k = 0, placed at the origin (z, z) = 0. (Since this is at zero 
momentum we cannot write this in the gauge (j3J).) The perturbative expansion is 
again expressed in terms of an ensemble of U(N) matrices with a finite convergence 
radius as before [9,12]. The series summation gives, for example, the A 00 component 



which can then be analytically continued to late times. 

In general, all closed string observables can be computed using the open string 
worldsheet theory, using similar perturbative expansions in terms of an ensemble of 
matrix models with a finite convergence radius, and the sum can be analytically 
continued to late times. For example, one-point functions of closed string on-shell 
vertex operators correspond to probabilities to produce closed string states in the 
decay [19]. The functions obtained can also be interpreted as coefficients of the 
associated component in the decaying brane boundary state [8,20]. The resulting 
closed string description is then valid at all times, even if the open string formalism 
is strictly valid only within the convergence radius, at suitably early times. 

2.2 The half S-brane as a Dyson Gas: partition function 

We can now show that all the S-brane disk amplitudes can also be expressed in terms 
of the correlation functions in a simple statistical mechanical system - the classical 
Coulomb gas of infinitely heavy unit charges on a unit circle, also called the Dyson 
gas. In particular, the S-brane disk partition function is related to the grand canonical 
partition function of the Dyson gas (see Appendix D of [12] for a brief discussion). 

To see this, recall first that a Dyson gas [21] consists of heavy unit charges confined 
to live on a unit circle in a two-dimensional plane. Pairs of charges, having positions 
e lti , interact via the two-dimensional Coulomb potential which is logarithmic. The 
kinetic term is zero since the particles are taken to be infinitely heavy and they do 
not move. The repulsive interactions give a two-body potential 




(15) 




(16) 



The grand canonical partition function of the Dyson gas is 



oo 



N 




(17) 



N=0 
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where z is the fugacity and Z^{j3) is the canonical partition function for N particles as 
a function of the inverse temperature (3 = 1/T. Putting in the logarithmic potential 
(ITtjj) . we see that at the special temperature /3' 1 = 1/2, 



00 7 N r N m 

N=0 J i=l 



n 

A<3 



ffti I 



(19) 



00 y Ar 



= EW rfc/ = r^' (20) 

where we used the identification of the integral over U(N) as in (jHJ) and (El). Thus we 
see that for (3 = 2, the interactions between the Dyson gas particles reproduce the 
integrations over U(N) in (jSJ). 

Comparing (ITOl) and ( 1201) it is tempting to identify z = —2n\e x ° and Z G = 
Z open . However, this is inconvenient because fugacity is the exponential of a chemical 
potential, and hence this identification would require an imaginary component for the 
chemical potential in our Dyson gas. Rather, Zq is more naturally identified with 
another quantity in the decaying brane worldsheet theory, namely the expectation 
value of an operator (—1) Nt , where Nt counts the number of tachyon vertex operators 
from the rolling tachyon background: 



1 



N 

— xj ^-*/i/\c j 

A --n 

(21) 



1 - 2n\e xU 



Identifying Z G «-> Z' open gives 



z = e^ = 2n\e x ° , (22) 

so the chemical potential \x essentially corresponds to target space time x°. Thus, con- 
sidering the decaying brane at different times corresponds to considering the Dyson 
gas at different points of thermodynamic equilibrium, corresponding to different val- 
ues of the chemical potential in the grand canonical ensemble, 

The interesting physical quantities in the target space of the decaying brane are 
computed as expectation values of closed string vertex operators with respect to the 
disk partition function Z open rather than Z' open . In order to obtain these observables 
from the Dyson gas formulation, we will need to compute the grand canonical partition 
function Z G as a function of a conventional fugacity and then analytically continue 

z -> ze tS -> ze tn = -z (23) 

to negative values of z, while continuing to identify time x° via the absolute value 
of z as in (1221) . Recall that because the late time physics (x° — > 00) lies beyond 
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the convergence radii of the series in (Q and (I20I) . an analytic continuation in the 
complex fugacity plane was in any case required. The continuation (I23p relating the 
string partition function Z open to the Dyson gas partition function Z G is simply an 
extension of this procedure. 



2.3 The half S-brane as a Dyson Gas: correlation functions 

Above we showed how the open string disk partition function is obtained from the 
thermodynamics of the Dyson gas. A general closed string correlation function is 
obtained by inserting vertex operators in the bulk of the disk. For homogenous brane 
decay, as described by ([1]), the spatial fields X\ i = 1, . . . , 25 in the vertex operator 
will simply contract amongst themselves and give a standard result on the disk. Thus, 
as in we will leave out this trivial spatial part of the computation since we are 
interested in the dependence of correlators on target space time. Then, the relevant 
part of the general vertex operator can be written as sums of products of terms of 
the general form 

: d ni X (z)d n2 X°(z)e ikoX ° : . (24) 

The correlation functions of such operators can be computed in the Dyson gas ap- 
proach by inserting external charges at the locations of the vertex operators and 
computing the expectation values of suitable moments of the potential. 

To illustrate, consider the stress tensor (n x = n 2 = 1 and k° — > 0). The T 00 
component contains the disk partition function Z open = B(x°) and the expectation 
value A 00 = 2(: dX°dX° :) ( fl2|) . In the Dyson gas thermodynamics, the latter 
quantity is computed by placing an additional charge k° at (z , z ) on the plane, 
inside the circle. The charged particles of the Dyson gas on the unit circle will 
produce an electrostatic potential acting on the bulk charge, 

^bulk-Dyson = _ fc ^ ^ |^ _ jt k |2 _ ^ 
k 

The holomorphic and anti-holomorphic derivatives of V"* 3 " 11 ^ ^ 011 will be called holo- 
morphic and antiholomorphic "forces" , 

a\/~huYk— Dyson »T^bulk— Dyson 
^ibulk— Dyson ^ibulk— Dyson (26) 

dz ' dz 

To compute ^4 00 in the Dyson gas, compute the grand canonical expectation value of 
the product of these "forces" , 

00 7 n r N Ht 

((F(zo)F(zo)))fi = Em n^^^t^ + ^H] 

n=o ' i=l 

. F bulk-Dyson (;Z0) {^^bulk-Dyson^ ^. }) ^ 
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then move the external charge to the origin, z = z = 0, and continue in the complex 
fugacity plane 

z -> -z = -27r\e x ° , (28) 
and replac^fl k° —> ik° to finally obtain 

((F(z )F(z )))p — > (: <9X°clX° :> = I • „4 00 (a; ) . (29) 

This is a prescription for computing the target space stress tensor for a decaying 
brane in terms of expectation values in the grand canonical ensemble for a Dyson 
gas. Here 

y Dyson = _ ^ ^ | _ g «, | ^ 

k<l 

is the potential energy of the Dyson gas. 

It is straightforward to generalize this to a comparison of more general thermal 
expectation values and more complicated closed and open string correlation functions. 
We add external charges at the locations of the vertex operators and compute the 
grand canonical expectation value of an appropriate moment of V. The potential 
between the external charges will produce an overall worldsheet position dependence 
in the amplitude (e.g., the initial factors in (J7])), while the interactions with the Dyson 
gas particles produce the part of the amplitude that is affected by the decay of the 
brane (e.g., the final factor F in (J7|)). Finally, to compute a string amplitude, the 
positions of all but one of the vertex operators are integrated over the diskfl This 
treatment can be applied to open string vertex operators too, except that here the 
external charges are inserted into the circle on which the Dyson gas particles are 
already present. The prescription is illustrated by some examples in the table below. 

3 The flow of time 

In the Dyson gas formulation, target space time has been replaced by the chemical 
potential \x of a statistical system. Different points of thermodynamic equilibrium 
characterized by different values of /i (keeping the temperature and the volume fixed) 
correspond to the different instants in time for the decaying brane. 

3.1 Time as the average particle number 

It is now simple to give an interpretation of time as a measure of the average particle 
number in the Dyson gas. In terms of the fugacity z, the average particle number is: 

N = zd z \nZ G = -^— . (31) 
x z 

3 Charges must be continued to imaginary values in order to make contact with real energies. 
4 We need 3 real parameters to fix the residual SL(2,R) symmetry. 



9 



In Rolling Tachyon CFT 


In Dyson Gas Ensemble 


. e i P °X°(T) . 

(on the boundary) 


Insert external particle with charge p° on the 
circle and modify the potential term by adding an 
interaction term for this external 

r*n T*crf>H t» q ri~ 1 r*l p qtiH T l^rcon era g 

L/lldl fci L.L1 Udl IjlUlC CI 11 LI LJ y oUll pLdO 


: d n X (r) : 


Insert charged particle on the 
boundary and compute the thermal average of 
the n-th moment of the potential in grand 

(•'.uimii i ■ - i t~i q y-j- Tn n f^i" inn r\j f hp 1 ni'cnn cr q g 
L^dllUlllL^dl JJdl 11 LlUll 1 HUll Ul tile l-'J'o'Jll g,cxo 


. e ik°X°(z,z) . 

(in the bulk) 


Insert particle with charge k° on 2d plane 
inside the circle at (z, z) and modify the potential 
by adding interaction between 
the bulk charge and Dyson gas particles 


: d ni X°(z)d n2 X°(z) : 


In cprf e> v"l~ f^i* Ti 1 r*n a ycrfiri t» i cA & 
illoCl Lj C-A.ljCllldl L/lldl tiCLl Udl tlL/1" 

inside the circle and compute thermal average of 
ni-th holomorphic and n2-th anti-holomorphic 

moments of the potential in the grand 
canonical partition function of the Dyson gas 



Table 1: Correspondence between state in CFT and thermal average in Dyson gas. 
Left hand side of the table is to be evaluated inside a correlator in rolling tachyon 
background. All the vertex operators in bosonic string theory can be obtained by 
taking suitable combinations of the above set. 



Recall now that every positive fugacity value was in one-to-one correspondence with 
an instant of time x° for the decaying brane, by 



\z\ = 2n\e x ° . (32) 

So the average number of particles in the Dyson gas related to the instant x° is given 
by 

N= 2?rA f =► s° = ln(^LVh27rA. (33) 

In particular, the infinite past x° = — oo corresponds to vanishing fugacity and N. The 
average particle number then increases monotonically as a function of the fugacity, 
corresponding to later time values x°. This relation is valid up to the point z — 1 
where N diverges. This corresponds to time value x° = — ln27rA. But we have also 
seen that this is also the boundary of the convergence radius for typical open string 
worldsheet calculations, beyond which it is more appropriate to use the closed string 
description of the system. 
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The average particle number N is a continuous quantity, just like time. Inter- 
estingly, however, the underlying physical quantity in the Dyson gas is the actual 
number of particles N. This is discrete, and fluctuates around N. The relative size 
of the fluctuations declines as N increases: 



At early times (small N), the relative fluctuations of A" are large and at later times 
N becomes more sharply defined. It is tempting to interpret this as "a continuous 
time emerging from an underlying discrete variable in the large A" limit" E 



3.2 A thermodynamic arrow of time 

Consider the grand potential of the Dyson gas ensemble, 

Q(fi,T,V) = -T\nZ G (fi,T,V) , (35) 

where we set Boltzmann's constant to 1. Q satisfies the thermodynamic identity 

Q — E — TS — fiN , (36) 

where E is the mean energy of the system. A Legendre transformation gives the 
conventional Helmholtz free energy: 

A(N,T,V) =n + nN = E-TS , (37) 

where we substitute for the chemical potential \x as a function of T, A" by inverting 
( }3~Tj) . A measures the free energy of the system as a function of the particle number 
(as opposed to the chemical potential). Setting T = 1/2 for the Dyson gas, we find 
that 

A(N, T,V) = — [(N + 1) ln(JV + 1) - iVln N] . (38) 

This resembles the familiar formula for entropy as a function of the mean number of 
particles. Indeed, the Shannon entropy in the particle number distribution is 

I s = - VVn) lnp(n) ; p(n) = — . (39) 

where p(n) is the probability that there are n charged particles in the Dyson gas. 
Doing the sum (1391) and using (13 ip for the mean particle number, we find 

I s = (N + 1) \n(N + 1) — Ann A" A(N, T, V) = - T I s ■ (40) 

The Shannon entropy increases with A^ and hence the free energy decreases with the 
mean particle number. Recall that the increase of A^ marks the passage of time. 
Thus the free energy decreases with the passage of time. We interpret this as a 
thermodynamic arrow of time for decaying brane universes. 

5 Quantitatively, one can cut off the infinite series summations at some value N max 3> N and get 
good approximations to all calculations. The size of the cutoff N max increases with time. 
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3.3 Equations of motion: generalities 

Usually, the flow of time is discussed in terms of a time evolution equation that 
determines the configuration of a system at all times, given an initial condition. 
Hence, it is interesting to ask how the spacetime equations of motion appear in the 
Dyson gas formulation of decaying branes. 

In string theory the spacetime fields appear as coupling constants in the worldsheet 
sigma model, and their equation of motion arises by requiring vanishing of all the beta 
functions [22]. In the presence of an unstable brane, the worldsheet theory includes 
a tachyonic boundary deformation. It is more subtle to derive the beta function in 
the presence of such a source, there are additional corrections which can be viewed to 
arise from higher order string diagrams in a degenerate limit [23,24]. The vanishing 
of the corresponding beta function then leads to a solution of the spacetime equations 
of motion containing a decaying brane. Specifically, the exactly marginal boundary 
deformation ([1]) preserves worldsheet conformal invariance and leads to a specific 
time dependent source term in the spacetime equations of motion. For example, the 
dilaton field equation in the presence of the decaying brane is 

V 2 = Jaa = (V^(0)e- ss ^y) disk , (41) 

where the r.h.s. is the one-point function of the dilaton vertex operator in the presence 
of the boundary deformation. For the homogenously rolling tachyon and a decaying 
p-brane in bosonic theory, the field equation turns out to be [25] 

V MN d M d N cl)(x , x) = c ■ 5 25 - p (x)f(x°) , (42) 

where M, N = 0, . . . , 25, c is a numerical constant and / is related to the target space 
stress tensor. Focusing on the simplest case of a space- filling brane (p = 25), the 
equation simplifies to 

-d 2 <j ) (x )=cf(x ) . (43) 

Since we have essentially reformulated worldsheet perturbation theory in terms of 
a statistical mechanical system, it should be possible to interpret equations like (1431) 
in terms of Dyson gas thermodynamics. Recalling that the time x° is related to the 
chemical potential as x° ~ 2/i, and that / is interpreted in terms of a thermodynamic 
average of moments of the electrostatic potential between the gas charges and an 
external charge at the origin, we expect an equation in the Dyson gas of the form 

- flj0(/i) = 4c • /(/i) • (44) 

In other words, if we identify the quantity <fr as an object in the Dyson gas, its value 
at all points of equilibrium (fixed by different values of p) can be determined via (j44D 
in terms of an "initial value" <f)(po) and the initial gradient d^<f)(po). 



12 



The identification of the spacetime fields within the Dyson gas formulation and 
the derivation of the equations of motion is done below and in the Appendix by 
translating the worldsheet beta function formulation into our statistical language. 
We have been unable to find a clear account in the literature of the derivation of 
the beta function equations for bulk fields in the presence of a source created by a 
tachyonic boundary deformation. While a complete calculation would clearly be an 
interesting basic string theory problem, the details are tangential to the discussion 
here. We outline the steps of the calculation in the Appendix. 



3.4 Equations of motion: derivation 

We begin by recalling how the field equations arise from requiring Weyl invariance of 
the decaying brane worldsheet sigma model. We start with the worldsheet action 

S = Sp + dSbdry 

= — *— f d 2 aVhh ab r] MN d a X M d b X N + A fds e x ° . (45) 
Aira' J J 

Next, we compute the disk partition function, as in (ITU1) . but now the contractions be- 
tween operators at the boundary will involve the curved space generalization dh(si, Sj) 
of the flat space distance function d[e %t \ e ltj ) = \e %ti — e ltj \. The correlators between 
bosons on the boundary become 

(X ( Si )X ( Sj )) = A^, Sj ) = ]nd 2 h (s h Sj ) . (46) 



The disk partition function then becomes 

-,x°\N 

J disk 

N=0 " ' ' " i ~" i<j 



Z.p x f /II,: II > (A^IV,!) , (47) 



where Z sp is the trivial contribution from the 25 free spacelike bosons X 1 . 

The corresponding statistical mechanical system (in addition to the free boson 
CFT) is the curved space generalization of the Dyson gas, a system where unit point 
charges on a (locally scaled) circle interact through a two-body Coulomb potential in 
curved space, 

^ Dyson = -^E A " dry ( s -^)' ( 48 ) 

i<j 

with the canonical partition function 

r N d, 

Z »W= /n|r ex P(-^ Dyson )- (49) 

•t i=l 71 

The inverse temperature is again (3 = 2. In the curved case there is no obvious ran- 
dom matrix interpretation. As before, we must consider a grand canonical ensemble 
partition function to make contact with the worldsheet theory. 
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The next step in deriving field equations is to introduce background fields in the 
worldsheet sigma model, for example a general target space metric Gmn(X) and a 
dilaton 4>(X), 

S[G, $] = — *— / d 2 a Vh [h ab G MN (X) d a X M d b X N + a'R$(X)\ , (50) 
4ira J 

where R is the worldsheet Ricci scalar. We take the target space metric to be almost 
flat, 

Gmn(X) = rjMN + £mn(X) (51) 
and choose €mn{X) to be a monochromatic gravitational plane wave 

e MN {X) = -^g c e MN e ik - x , (52) 

where g c denotes the string coupling. Similarly, we take a plane wave of the dilaton: 

$(X) = -A*Kg c <Pe ik - x . (53) 

Next, perform a Taylor expansion in the metric and dilaton perturbation for the disk 
partition function, so that it becomes 

Zdisk[£, 0] = (e- SSbd ^ {l + g c : V : +|| : VV : + ■ • • }>„ (54) 

where 

V = - [ d 2 aJh{a) \e M N h ab (a) d a X M d b X N e lk ' x + a'tRe*'*] (55) 
a' J I J 

is the graviton-dilaton vertex operator. Requiring Zdisk to be invariant under the 
Weyl transformations h a b — > h' ab = e 2w ^'h a b then imposes the on-shell conditions 
and physical polarization conditions for e MJV . Or, more precisely, in coordinate space 
the /?-function equation gives the (linearized) Einstein equation (see e.g. [26]) in the 
presence of the rolling tachyon source in string framed 

R L MN + 2V L M V L N cj ) = ^Tt IN . (56) 

Here the superscript L means linearized. Recall that the gravitational field from 
which the Ricci tensor is constructed is (IBTj) ; to first order in e 

r mn = \ {d 2 t MN - d P d M e p N - d P d N e p M + d M d N e p P } . (57) 



Meanwhile, the stress tensor Tmn appearing on the right side of (15611 is equal to the 
disk one-point functions (fl2"]) - (TT3l . 



3 The derivation of the source term is subtle. See the Appendix. 
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We now have sufficient information to identify what computation in the Dyson 
gas gives the thermodynamic counterpart of the Einstein equation (1561) . Generating 
a weakly curved background spacetime metric corresponds to introducing additional 
bulk charges to the system, and considering the thermal expectation values of some 
(generalized) forces resulting from the point charges on the circle acting on the added 
test charge. Then, just like above, we should require scale invariance to derive the 
counterpart of the field equations. It is possible that the requirement of scale invari- 
ance is tantamount to requiring that the statistical system be at a critical point, but 
we have not explored this carefully. In any case, deformations of the Dyson gas which 
introduce additional external charges result in deformations of the spacetime. 

In order to keep the exposition brief, we will ignore the 25 free spacelike bosons 
that play role in determining the full Einstein equation (j56p . Thus, in what follows 
we will only focus on the time-time component 

R 00 + 2d 2 Q <P = 8vrT 00 , (58) 

as it will be straightforward to repeat the calculations for the other components. 
Our starting point is fl54|) . which instructs us to consider the quantity 



Z« disk »M] = Z sp x (((i)) p + gc ((0)) fi + ...) , (59) 

where Z sp denotes the contribution from the 25 spatial directions of (j5"%|) . and the rest 

is the non-trivial contribution from the grand canonical ensemble of charges. We can 

calculate that the second term in the expansion, corresponding to the second term 
( e -ss hdiy j dk v ^ in (|5^ is 

« a )^ E/ EM ((aw(fc0)) ^' (60) 

^=0 



where 



((d N (k°)))^ N ee -e 00 / d 2 a^h(a) / J] p exp { - ft [ : ^ bulk (a,a') : + 



x 



V -Dyson . T rbulk— Dyson 
h + V h J 



+ 



J uabf) T/ bulk-Dyson o T/ bulk-Dyson J- iab . a a T /bulk . I 

-—h d a V h d b V h — — — n .d a d b V h . | CT ^ CT 



(A; ) 2 ° 11 ° h 2k° 

In the above, o = (cr 1 , a 2 ) is the position of the bulk charge inside the disk, e ab is the 
unique antisymmetric tensor in two dimensions and 

. T^bulk/ /\ I i T rbulk— Dyson 
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is a regularized^ net potential energy measured at a due to a bulk charge k° placed 
at er and Dyson gas on the boundary of the unit disk, subject to the Neumann 
boundary condition on the boundary. These terms can be expressed in terms of curved 
space generalization of bulk-bulk (X°(a)X°(a')) and bulk-boundary (X°(a)X°(si)) 
potentials respectively: 

^bulk-Dyson = _ ^0 ^ A bulk-bdry ^ ^ 

i 

VTW) = -k°A™ k (a,a'), (61) 

where both A^ ulk-bdry (<7, Si) and A^ ulk (cr, a') are subject to the Neumann condition 
on the boundary. The expression inside the square bracket of flSTj) is the product of 
the generalized forces on the bulk charge from the Dyson gas. Those expressions are 
regularized since these are measured at the location of the bulk charge itself. We 
denote this regularization scheme by : :. Finally, to get the expression of (1541 from 
(l6Tj) . we follow these two rule^f): (1) Analytically continue z — > — z = —2n\e x °, (2) 
Replace k° by ik°. 

In the thermodynamical system, the above quantities are thermal expectation 
values of particular moments of the Coulomb interaction potentials. Requiring 0591) 
to be scale invariant yields the Einstein equation fl56|) . (For additional details that will 
be involved in the derivation, see the Appendix.) Derived in this way, the Einstein 
equation does not refer to spacetime tensors; rather it refers to Fourier transformations 
of quantities in the momentum space (k°,k). This is related to the familiar Ricci 
tensor by the following steps. First, Fourier transform the Ricci tensor of 036]) from 
coordinate space to momentum space: 



RMN(k) = J dx° J d 25 x e ik ' x R M N(x) 



I 

= -{k 2 e M N(k) - k P k M e p N - k P k N e p M + k M k N e p P } . (62) 

Then, do the inverse Fourier transformation, but now with the chemical potential /i 
instead of time x°. In particular, this defines 

R m (fi,x) =J^J d 25 ke- tk °^R 00 (k ,k) . (63) 



In our example, the time-time component of the Einstein equation 0581) is con- 
verted to an equation for thermal expectation values 

+ 2dl<P = 8nT^ , (64) 



7 We subtract the usual divergent selfinteraction contribution of the bulk charge but keep the 
finite contribution from its image charge necessitated by the Neumann boundary condition. 

8 Note also that the spectator CFT contained in Z sp is required to fully recover the answer of 

mi. 
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where fj, = In z/f3 is the chemical potential, and the stress tensor component T MAt is 
the thermal expectation value (1271) in the limit k° — > 0. This particular component 
of the stress tensor turns out to be a constant, the conserved initial energy. Non- 
trivial chemical potential dependence would be obtained for what where initially the 
pressure components T„ of the stress tensor. 

We have sketched how the thermal ensemble counterparts of the field equations 
emerge from requiring scale invariance of a Dyson gas with additional test charge 
density included. In this way we could also derive equations of motion of other 
spacetime fields, such as the dilaton equation which we discussed in section 3.3. 

4 Discussion 

In this paper we reformulated the worldsheet description of brane decay in terms of 
the statistical mechanics of the Dyson gas. The progress of time was marked by an 
increase in the average number of particles in the gas. This gave rise to a thermody- 
namic arrow of time - the decrease of free energy (or, equivalently, the increase of the 
Shannon entropy of the particle number distribution) marked the passage of time. 
The spacetime equations of motion translated into relations between the expectation 
values of different statistical moments at different values of the chemical potential. 
Our analysis has essentially involved rewriting worldsheet perturbation theory in a 
statistical language^, and we are not precisely describing the "emergence" of time 
from a timeless system. Nevertheless, there are several lessons to be learned here 
concerning what it might mean for time to be emergent. 

First, in any scenario where time is emergent, the absence of time in the under- 
lying microscopic description will require the latter to be some kind of Euclidean or 
statistical system. In analogy with the AdS/CFT correspondence one might have 
expected some kind of relation between larger scales in the statistical system and 
later instants of time. In the system we have studied, the "scale" in question turned 
out to be the number of interacting particles in the system, and time as a continuous 
variable arose from this discrete quantity via a Legendre transform relating it to a 
chemical potential. 

Second, in a scenario where time emerges from a statistical system, one might have 
expected a non-equilibrium flow to realize the flow in time. In fact, this is a somewhat 
misleading expectation. In non-equilibrium statistical mechanics, there is already a 
physical time and the dynamical equations describe the flow in time of statistical 
moments towards their equilibrium values. In our setup, we instead found that each 
instant of time was described by a different point of equilibrium in a statistical system. 

9 The electrostatics interpretation of worldsheet Green functions as 2d Coulomb potentials be- 
tween point charges is straightforward. The main issue was that here it specifically leads to a 
statistical mechanical system with a reinterpretation for time in target space. 



17 



The spacetime equations of motion became relations between these different points 
of equilibrium. Indeed, it is possible that any emergent or holographic description of 
time will have such a character. We are attached to the idea of a "flow in time" partly 
because it is conventional to construct the initial value problem in terms of data on a 
spacelike surface. But we could equally well have formulated it on a timelike surface 
of co-dimension one. Indeed, the latter is the natural way to do things in AdS space. 

Most of this paper was concerned with a statistical formulation of the rolling 
tachyon background of Sen ([1]). However, in order to derive the equations of motion 
we required Weyl invariance of deformations of the Dyson gas. These deformations, 
involving the addition of external charges, gave rise to new spacetime backgrounds 
solving the equations of motion. Thus the set of Weyl invariant deformations of the 
Dyson gas with (3 = 2 all give statistical descriptions of spacetimes with decaying 
branes and different initial conditions^]. This is reminiscent of the AdS/CFT corre- 
spondence where deformations of the CFT are dual to deformations of AdS spacetime. 
From this perspective, one difference is that in AdS/CFT the deformations need not 
be marginal - relevant deformations are dual to asymptotically AdS spacetimes that 
differ in the deep interior, and the RG flow of the field theory is related to the space- 
time equations of motion. 

A simple way of exploring the role of non-marginal deformations of the Dyson gas 
is to consider generalizations of Sen's background (pQ) that are of the form 

SS open = \Jdt e aX °® . (65) 

The parameter a in the exponent is related to the inverse temperature of the associ- 
ated Dyson gas (3 by a = ■sj (3/2. For j3 < 2, the deformation is relevant, for (3 > 2 it 
is irrelevant. This general system has been studied in the random matrix literature, 
and is known as the circular /3-ensembli"1. Dyson conjectured that the canonical 
partition function of the ensemble at fixed iV is 

Z N ((3) = y ) . (66) 

K ^ J [r(i + f)p 1 ; 

and various proofs have been presented in the literature (see [28]). Using Dyson's 
formula, we can construct the grand canonical partition function 

oo 

Z g (z,(3,V) = J24 Z n(P)/N\ , (67) 

N=0 



10 Note that although we considered deformations at the boundary of the worldsheet correspond- 
ing to the one-dimensional Dyson gas, this is just a specific example of a more general relation 
between scale invariant statistical mechanical systems and (time-dependent) string backgrounds. 
For example, it would be interesting to find a generalization for closed string tachyon condensation. 

11 Incidentally, a longstanding problem was to find a random matrix model formulation of the 
general /3-ensembles. For circular ensembles, a solution in terms of matrix models of certain sparse 
unitary matrices was found in [27]. 
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with Z N (j3) given by (1661) . It is readily shown that the partition function diverges 
for (3 > 2 and converges for (3 < 2, for all values of fugacity z. Exactly at (5 — 2, 
the sum is convergent for \z\ < 1, and requires analytic continuation to larger values 
of z as described before. This suggests that (3 = 2 is a critical point for the Dyson 
gas (see [29,30]). Furthermore, the renormalization group flow [31] of the relevant 
deformations (j3 < 2, convergent partition function) could perhaps be related to the 
spacetime equations of motion as in AdS/CFT and the worldsheet analyses of [32]. 
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A Toward deriving equations of motion using Weyl 
variation method 

A.l Goal 

We would like to show that the equation of motions of the bulk massless fields, like 
graviton, dilaton etc. in presence of a decaying brane can be derived using Weyl 
invariance of the worldsheet non-linear sigma model. In absence of any boundary 
perturbation, the derivation is easy and given in Polchinski's book [26]. In presence 
of a boundary perturbation, the derivation is not so straightforward. 

The sigma model is described on a curved worldsheet in presence of massless 
closed string fields background and boundary perturbation corresponding to decaying 
brane. In particular, for the graviton-dilaton we would like to show that its equation 
of motion is given by 

R L MN + 2V L M V L N cj> = T MN , (68) 

where Rmn ls the linearized (in Iimn) target spacetime Ricci tensor and Tmn is the 
graviton-dilaton one-point function in presence of boundary perturbation , dS^y = 
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\fdse x "^ : 

T M M = (Ve-"f d - X% \ , (69) 

where V is the graviton-dilaton vertex operator defined in (I55p and (0) v = J T>XT>he~ 
implies that the correlation function is evaluated in the free theory using perturbative 
methods. Here Sp is the Polyakov action defined in the first term of (HoT) . 



A. 2 Steps 

How can we derive (1681) ? By enforcing Weyl invariance in the full quantum theory of 
the worldsheet non-linear sigma model, we can achieve this. In this method we shall 
be able to derive the analog of fl68|) in momentum space. We give the steps to be 
carried out in this section : 

1. The starting point is We write Z disk [e, <fi] = Z + Zi + . . ., where Z n is the 
term with n-th power of V : Z n = ^ e ~' 55bdry (: V :)"%• Since SS^dvy involves 
exactly marginal boundary perturbation, Z is invariant under worldsheet Weyl 
transformation : 5\yh a b = 5u(a)h a b. However, other terms Z n {n > 0) are not. 
If we make sure that Z\ is Weyl invariant, all other terms with higher powers 
of V will automatically be so. The Weyl variation of Z\ is (the Polyakov action 
is Weyl invariant): 



b w Z x = J VXVh e- Sp 5 w : V : 



+ lvXVhe^J25 w ^- X $ dS N f (S):)N :V:) . (70) 



N=l 



Weyl invariance is imposed by demanding 5y/Zi = 0. The notation : : around 
an operator implies that they are normal ordered. 

2. The first term in the r.h.s. of the above equation can be derived by using 
the result given in Polchinski's book. On curved worldsheet, first, we define 
normalized bulk operator like V as: 

: V := e^V, D Mk = / fafaAn.,.) ^^^^ . (7!) 

Here A^ ulk (o"i, <r 2 ) is the bulk-bulk propagator on curved worldsheet. Note that 
in the flat worldsheet limit (+ boundary conditions) 

A£ u Vi, <t 3 ) - A^, z; w, w) = \ In \z - w\ 2 . (72) 
Next we define the Weyl variation of the normalized bulk operator as 

5 W : V := e SwBb ^V + 6 C ^ P V , (73) 



20 



where the last term in the above equation denotes the explicit Weyl variation. 
The Weyl variation of D bu i k is 

^D bulk = / dVA^A^Vi,^)-— S ——±— . (74) 

3. Choosing a particular RG schemed on the worldsheet, it turns out that 
6 w :V:=g c J d 2 aVhSco R'£ IN (k)h ab : d a X M d b X N e ik ' x : +a'FR : e 



JkX 



(75) 



where 



R'mnW = ^{-k 2 e M N + k N k p eMP + k M k p s N p-k M k N e P p + Ak M k N (f)} 
= RM N {k) + 2k M k N (f) 
F = -\k 2 ^. (76) 



4. The most non-trivial part for evaluating the r.h.s. of fjTOl) is the second term. 

We need to generalize the rule in (1731 to the case when we have a composite 
operator made of bulk and boundary operators on the curved worldsheet. We 
would like to give a recipe for this generalized definition of normal ordering 
next. 

5. Polchinski's Rule Generalized : We use <Tj = (of, of,) for bulk and (t^Uj) 
for boundary coordinates. We define the following two operators Dbdry and 
Dbulk-bdry on the boundary and bulk-boundary respectively, 

Dbdry EE \ J d Ul du 2 A bdry fo, U 2 ) 

1 f s s 

Dbulk-bdry = o / d (T 3 du 3 A bulk _ bdry (cr 3 , M 3 ) — - — - . (77) 

2 J 0A^(cr 3 ) ()A p (M 3 ) 

In our case the normal ordered boundary operator is : e^ ^ := e E>bdry e x °^- ) 
which is exactly marginal. Hence its explicit Weyl variation as well as the Weyl 
variation of D bc j ry are vanishing. Keeping these two points in mind, we define 
the following generalized differential operator which does not include D b d ry : 

D = D bu i k + D bu i k _ bdry . (78) 

Let T be a composite operator made of bulk and boundary operators. In our 
case T will be A)e x ° (Sl - ) V. We define the Weyl variation of the normal 

ordered T as 

5 W : T := (D5 W D)*F . (79) 
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We put 7 = in eqns (3.6.17a-c) of [26]. 
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The subscript * denotes an important rule : we set terms like (D bu i k _ bdry ) n = 
for n > 1. The reason is that there is no concept of self-contraction for 
computing correlators between bulk and boundary operators. The explicit Weyl 
variation of T in our case is zero since the variation for both \\ i e x °( s i) anc i y 
are zero. Equation (I79|) will be our working definition. 



A. 2.1 Explicit calculations 

From (JZZ5 and (JZSD, we find that 

(D5iyD)*jF = (Dbulk<WDbulk + Dbulk<5tyDbiilk-bdry + Dbulk-bdry^iyDbulk) J~ ■ 

In our case 

T= f d 2 aVh\e M Nh ab d a X M d b X N e lkX (a) + a'R(Pe lk - x ] / ' JJ rfs fc JJ(-A)e 



SO) 



(81) 

For T given in (IHTj) . we compute its Weyl variation 5w '■ J 7 ', following (!79|) . So we 
need to evaluate (IBUl) . The result is : 

(D5 W D)* T 



d 2 a\\ds k Vhh ab e 



£MN 



- Ak M k N (<9 a (WAbuik<9bAbuik + <9 a A bulk <9 6 <VA bu ik) 



2fc 2 |r / MJV (9 a ^^Ab u i k )A bu i k + V MN (d a d' b A hnlk )5 w A hnlk 
+ ik^ M (d {a 8 w A hvil ^A hvilk d b) X N ^ 

-bdry(c", bulk 

i 

+ 2ik°rj MN 22 5H/A bu ik-bdry(o-, Si)d a & b A hxilk 

i 

+ 2tk° v MN d a d b 5 w A hulk J2^ui^cT,s l ) 



+ 2 i k iM r ] m) d {a 6 w A hulk Y,d b) A hnlk -bdry(c, Sj 



12) 



This equation actually evaluates Weyl variation of (Ve~ <5,Sbdry )/ l where the subscript 
h denotes curved worldsheet. This is easy to verify : first, we compute one-point 
function of V in presence of boundary perturbation on flat worldsheet, (Ve~' 55bdr y) rr 
Next, we covariantize this result for curved worldsheet to obtain (Ve~ SSbdiy )h- Finally 
we can compute its Weyl variation by using the rules in (!73|) . (ITBl) and (1791) . 

The next step is to express the r.h.s. of (182]) as ~ 5uT A l I . Finally we have to 
substitute (182|) . (!75|) and (!76|) in r.h.s. of (ITU]) to obtain the momentum space version 
of equation of motion in (IqT)1) . We leave this issue for future work. 
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